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Abstract 

We introduce a new approach to the Josephson effect in SIS tunnel junctions. 
The Josephson coupling energy is calculated from the overlap of real space 
Cooper pair wavefunctions in two superconductors through an insulating bar- 
rier. It is shown that the Josephson tunneling is limited by the size of the 
Cooper pair and its shrinking during the tunneling. Therefore, the Josephson 
coupling energy and the critical current become extremely small in high T c 
superconductors, including MgBi- This shrinking also causes the observed 
DC supercurrent in low T c superconductors, such as Pb and Sn, to fall off 
much faster than 1 /R n for tunneling resistance R n above several ohms. Con- 
sequently there is a material-dependent threshold resistance, above which the 
supercurrent decreases much faster with increasing resistance. The impurity- 
induced shrinking is also shown to limit the critical current. Furthermore, the 
(weak) temperature dependence of the Cooper pair size is found to contribute 
to the temperature dependence of the DC supercurrent. 
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I. INTRODUCTION 



The Josephson effect may be the most fascinating and intriguing property of 
superconductors. 1,2 In 1962 Josephson predicted that supercurrents can flow through the 
insulating barrier between two superconductors due to the Cooper pair tunneling. 1 The 
supercurrent depends on the relative phase of two superconductors, where the phase is as- 
sociated with the order parameter of the superconducting condensed state. The prediction 
was confirmed in experiments quickly 3,4 It is interesting that the experiment by Nicol et 
al. 5 showed already the possible DC supercurrent, although its reality was not noticed. This 
discovery led to the fabrication of the Josephson junction devices and the high precision mea- 
surement of fundamental constants. 6 Now it seems that our understanding of the Josephson 
effect is on firm ground. 6,7,8 

Josephson 1,2 used the tunneling Hamiltonian 9 to calculate tunneling currents between 
two superconductors and noticed that even at zero applied voltages, a DC supercurrent can 
occur. At finite voltages V, he found an AC supercurrent of the same amplitude as that 
of the DC supercurrent, with frequency 2eV/h. The amplitude of the supercurrent was 
calculated by Anderson 7 and Ambegaokar and Baratoff. 10 At T = OK their result of the DC 
supercurrent, j, is 

j = j lS iri(f), (1) 

with 

* A 

Here is the phase difference of the two superconductors with the same energy gap A. 
Note that the amplitude of the supercurrent is proportional to the energy gap and inversely 
proportional to the tunneling resistance. 

However, there are still some fundamental experiments which remain puzzling. 11,12,13,14,15 
For instance, the pair-quasi particle interference term was shown to be negative in experi- 
ments, whereas the theory predicts positive sign. 11 For AC Josephson effect, Likharev even 



claimed that there are more questions than answers. 12 Furthermore, high T c cuprate Joseph- 
son junctions are mainly SNS type junctions, since SIS junctions do not show the Josephson 
effect. 13 This result is not consistent with the above theory, which predicts bigger Joseph- 
son current for higher T c . 1,2 ' 7,8 ' 10 Recently, this behavior has also been found in MgB 2 SIS 
junctions. 14 ' 15 In particular, the Josephson supercurrent was not observed for the big gap 
(~ 7meV) of MgB 2 . We stress that this behavior is closely related to the anomalous de- 
pendence of the maximum DC Josephson current on the tunneling resistance R n , found in 
low T c superconductors, such as Pb and Sn: 16 ' 17 ' 18 ' 19 the Josephson current decreases much 
faster than l/R n above several Q. 

In this paper we present Cooper pair wavefunction approach to the Josephson effect. It 
is shown that the Josephson coupling energy is determined by the overlap of the Cooper 
pair wavefunctions of two superconductors divided by a thin insulating layer. (In fact, this 
idea was suggested by Josephson in his original paper, 1 and confirmed by Ambegaokar and 
Baratoff. 10 However, it was not pursued thoroughly.) The critical current is, then, calculated 
from the coupling energy Ej in the usual way. 2 Since the Cooper pair tunneling can be 
understood more easily from the Cooper pair wavefunction, this method is advantageous 
to deal with the Josephson effect in high T c superconductors including MgB 2 and for the 
insulating barriers with high tunneling resistance R n . We have found that the Cooper pair 
tunneling is strongly limited by the size of the Cooper pair and its reduction during the 
tunneling. Since the Cooper pair size is also reduced by the impurity potential scattering, 
ordinary impurity also limits the supercurrent. As a result, the Josephson current is much 
smaller than expected in higher T c SIS junctions and falls off much faster than 1/R n for 
thicker and/or higher insulating barriers, in good agreement with the experimental findings. 
The preliminary result was reported before. 20 
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II. COOPER PAIR WAVEFUNCTION APPROACH TO JOSEPHSON 

TUNNELING : T=OK 



A. Josephson Coupling Energy 

The approximate expression of the supercurrent suggested by Josephson 1 is 

3 = \jl^r + \hWh (3) 

where ipi and ip r are the effective superconducting wavefunctions on the left and right sides, 
respectively 21,22 Actually, it can be exact and it is better to write first the Josephson coupling 
energy Ej in terms of the Cooper pair wavefunctions and then calculate the supercurrent j 
from Ey? 

3 = "A (4) 

J h d<f) y J 

Figure 1 shows the real space Cooper pair wavefunctions of two superconductors with 
the same energy gap near the insulating barrier. The solid lines denote the Cooper pair 
wavefunctions in the absence of the barrier, whereas the thick lines denote the change of the 
Cooper pair wavefunctions due to the insulating barrier. The Josephson coupling energy 
is, therefore, determined by the overlap of the Cooper pair wavefunctions, connected by 
the phonon Green's function. 23,24 Previously, this idea was employed in some cases. 25,26 
For instance, for a pure superconductor without a barrier with the Einstein phonons, the 
coupling energy or the pairing energy is given by 

= V S 2 U kVkUk>V k > , (5) 

k k' 

where the (effective) Cooper pair wavefunctions F(x, y)(— Y,k Uk v k elk '^ r_y ^) and F*(x, y) are 
coupled by the Einstein phonon Green's function, i.e., the Dirac delta function, S(x — y). 
Here V is the phonon- mediated matrix element. In the same way, the Josephson coupling 
energy Ej is given by 



Ej = vJ J F r *(x,y)F l (x,y)5(x-y)dxdy + V J J F l *(x,y)F r (x,y)5(x - y)dxdy 

= V J F;(x)Fi(x)dx + V J F;(x)F r (x)dx, (6) 

with Fi and F r being the effective Cooper pair wavefunctions in left and right sides. Observe 
that this expression has the same form as the approximate one, Eq. (3) except the integra- 
tion. In fact, Ambegaokar and Baratoff 10 noted that Josephson coupling (free) energy can 
be calculated by the overlap of Cooper pair wavefunctions, which agrees with our result: 

Ej = f\d\/X) < \H T >, (7) 
Jo 

where A is an explicit coupling constant and Ht denotes the tunneling Hamiltonian. Figure 
1 shows that it is essential to calculate the tail of the Cooper pair wavefunctions to determine 
the Josephson coupling energy. 

It is necessary to emphasize that Eq. (6) leads to the familiar result in k-space. From 
the tunneling Hamiltonian 9 

H T = ^2(T kk ,cp- s c l ks + H.C.) (8) 



kk's 



one finds the Josephson pair tunneling Hamiltonian, Hj 27 ' 28 ' 29 



Hj = - E(^4KV-Xt + H - c -)- ( 9 ) 

kq 

Then, the coupling energy is 

e.j = - j(E u >k E u \ v \ + E u >l E «) 

k Q k q 

—^(E^ + E^) 

kq kq 

with Jk q = J. Comparing Eqs. (6) and (10), it is evident that J is proportional to the 
phonon-mediated matrix element V, i.e., 

JocV. (11) 

Accordingly, it is crucial to note that J kq is the pair scattering matrix element across the 
barrier due to both the tunneling and the electron-phonon interaction. If we define 



^ r EE Y, F l = |^r|e i0r 
k 

^=E^ = l^|e^ (12) 
q 

Eq. (10) is rewritten 

Ej = -J(^% + = -2J|tt,||tt r |cos(& - r ). (13) 

It should be noticed that the result is formally the same as the internal Josephson 
coupling energy in two-band superconductors, 27 ' 28 and believed to reproduce the Josephson 
coupling energy Ej in conventional SIS junctions. 27 ' 28 ' 29 However, the exact expression of 
J kq ee J is missing, although the (2nd order) tunneling-induced pair transfer term (J/Q, 



jT = l Tfcgl (U) 

was confused with J. 27 ' 29 Here E l k and E r are quasi-particle energies. It would be dangerous 
indeed to suppose that J kq = J (phonon-mediated pair-scattering matrix element across the 
barrier) and J kq (tunneling-induced only pair-scattering matrix element) are the same when 
they are not. Note that Jj g is not proportional to V, unlike J, so Eq. (10) or Eq. (13) with 
J kq does not reproduce the previous result, Eq. (2), which remains puzzling. 29 
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B. k-space Approach 

Now we need to compute Eq. (6) or Eq. (10) to find the Josephson coupling energy Ej 
and the Josephson current j. It is clear that both equations should lead to the same result. 
First, we consider Eq. (10), which requires the calculation of the pair-scattering matrix 
element J kq . For simplicity, we assume two superconductors with the same gap, A, which 
will be generalized later. From Eq. (9) we obtain 

Jkq =< k t-k I \V e - ph \ -q [,q T>' (15) 

where V e - V h is the phonon-mediated electron-electron interaction. Tunneling effect is in- 
cluded in the states of the left and right sides. It is noteworthy that this situation is much 
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the same as the calculation of the matrix element between scattered-state pairs in Ander- 
son's theory of dirty superconductors, 30,31 where impurity effect is included in the scattered 
states and the electron-phonon interaction couples the states. 

Using the Tunneling Hamiltonian, Eq. (8), the one-particle wavefunctions are given by 

^(r) = JV s W(r) + E7^V«(0] 
k e i tk 

V k {r) = NMltr) + £ -^-#(r)] (16) 

where (fi l q and <j) r k are the left- and right-hand states 32 with e q and being the electron 
energies. N k and N q are the normalization constants. Prange 32 proposed the nonorthogonal 
almost confined states for the left- and right-hand states, whereas Bardeen 33 introduced the 
WKB wavefunction, 



(f) l m = Cp- 1/2 e i{p y y+PzZ) sin{ Px x + 7), x < x a (17) 
C\p x \- 1/2 e l(p y y+p * z) exp{- f X \p x \dx), x a < x < x b (18) 

Jx a 



m 2 



where C is a normalization constant and \p x \ = (2fj,U — p 2 y — p 2 ,) 1 ^ 2 with U(x) being the 
potential energy, u, is the electron mass. The barrier extends from x a to x b - Notice that 
the WKB wavefunction is a sinusoidal stationary wave (in x-direction). Thus, we have 
degenerate pairs m = (p x ,p y ,p z ) and fh = (p x , -p 2 ). 

For the Einstein phonon model, 23 ' 24 ' 31 the pair-scattering matrix element is 



J kq = v 



J dr^(rW-l(rW. q (rW q (r) 
= vjdr\r k (r)\ 2 \4(r)\ 2 , (19) 

which denotes the (density) correlation function between the eigenstates ip l q (r) and i^l(r). 
Since the Bardeen's wavefunction is not much different from a plane wave, we may employ 
plane waves for <p l q and <p r k to obtain the matrix element, 

(20) 
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where Qi and Q r are the volumes of the left and right sides. Unfortunately, the sums are 
divergent, which are common in the perturbation theory for continuous spectra. The remedy 
is to use the scattering theory with a cutoff (£ ), consistent with the bound state of Cooper 
pairs. 34 Since each sum denotes the relative probability contained in the virtual scattered 
wavelets (due to the barrier) within the BCS coherence length £ — hvp/nA (compared to 
the plane wave part), 34 we find 

1 V l T ^l 2 ~ c IT \*> Nf ^P -<\T \*> NF7r 



1 V l T ^l 2 ~ c IT \*> Nf7T % -<\T \*> Nf7 



and 



< J kq >= J « V < \T\ 2 > (22) 

Here N F is the density of states at the Fermi level and the angular brackets indicate the 
average value over the states. 

Consequently, substituting Eq. (22) to Eq. (13) the Josephson coupling energy is written 

2n 

Ej^-—< \T\ 2 > N 2 F A cosifa - <f> r ), (23) 
where A = N F V . Accordingly, the Josephson supercurrent is 

3 = t-^- sin((f>i - <f> r ), (24) 

with 

h = IA- (28) 

Although this equation for the maximum DC supercurrent is similar to the previous one, Eq. 
(2), physically it is significantly different. The factor 1/A shows that the Josephson coupling 
energy and the supercurrent depend on the superconductor. (Notice that this factor also 
appears in the initial T c decreases due to magnetic impurities 34 and weak localization. 31 ) 
As a result, if we include the reduction of the Cooper pair size during the tunneling, which 
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is important especially for superconductors with high T c , the resulting supercurrent is not 
similar at all. 

When the superconductors have different energy gaps, Ai and A 2 , the pair-scattering 
matrix elements are 

< J kq J « Vl < |T| 2 > ^ + V 2 < |T| 2 > ^ (26) 

where Np,i and N F:2 are the density of states of left and right sides, respectively. Subse- 
quently, one finds 

Ej = -7T < |T| 2 > iV ljF iV Fi2 A 1 A 2 (-^— + — L- )cos(0, - r ), (27) 

A1A1 A 2 IA 2 

and 

j = ^-e< |T| 2 > N FA N F;2 A 1 A 2 (-^— + - r ), (28) 

with 

i^hh^xk + ^J- (29) 

Here Ai and A2 are the BCS coupling constants for two superconductors. 



C. Cooper Pair Wavefunction Approach 

We consider now Eq. (6), the Josephson coupling energy in terms of the Cooper pair 
wavef unctions. This method is more appropriate for high T c superconductors, including 
MgB 2 , which have tightly-bound Cooper pairs with small size. It is also desirable for the 
insulating barrier with high tunneling resistance R n , since the Cooper pair size shrinks 
significantly during the tunneling in that case. A similar reduction of the Cooper pair size 
caused by the ordinary impurity scattering requires this approach, too. 

The Cooper pair wavefunctions are expressed in terms of the one-particle wavef unctions, 
Eq. (16): 
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F l (x)=J2u q v^ l q (x)^_ q (x) (30) 
q 

F r(x) = J2 U kVkVk( x W-k( X )- ( 31 ) 

k 

Substituting Eqs. (30) and (31) into Eq. (6), it is straightforward to show that 

E > = (S? + 5? dS?'?^ + (32) 

which is indeed the same as that obtained by the k-space approach, i.e., Eq. (10) with Eq. 
(20). Whereas this real space approach allows us to take into account the change of the 
Cooper pair size during the tunneling through the Cooper pair wavef unctions. 

If we insert the Bardeen's WKB wavefunction into Eq. (30), the left-side Cooper pair 
wavefunction shows the exponential decay in the barrier region: 

Fi(x) = ^C£u p v p )\p° x \- 1 e- 2 ^-^ x a <x<x b (33) 
v 

where = hn = (2/j,U — p 2 — pi) 1 / 2 with U being the potential energy. Beyond x b , the 
tail of the Cooper pair wavefunction should be, after utilizing the one-particle wavefunction, 
Eq. (16), 



\T, 
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q k \ 6 q ek > 

We have used the fact that only (<f) r k , (ff_k) pair contributes to the coupling energy, Eq. (32). 
(It is obvious that \T kq \ 2 oc e~ 2Kd , with the barrier thickness d = x b — x a . 35 ' 36 ) However, since 
(f) r k and 4> r -k are basically plane- wave-like, Eq. (34) shows not the presumed exponential tail 
but a constant amplitude. Of course, the divergence in the sum implies this behavior needs 
to be corrected according to the bounded Cooper pair wavefunction. Therefore, we assume 

F l (x)^<\T kq \ 2 >J2u q v q ^f^e-^- x ^ x>x b . (35) 

(In fact, x denotes the center of mass coordinate of the Cooper pair rather than the relative 
coordinate. Since the phonon-mediated attraction between the electrons in the left-hand 
side can propagate only up to the range of the Cooper pair size £ m the right-hand side, 
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this assumption seems to be reasonable.) Inserting Eq. (35) into the first integral in Eq. (6) 
we obtain (in the right side) 

V [ r F r *(x)F l (x)dx = V<\T\ 2 > ]>>X*E<^^ 1^ e ~ {X ~ Xb),i ° dx 

J k q aiIVf J Xb 

= y<iT| 2 >E^rE«^ (36) 

k g 

where J r denotes the integration over the right side and x a = 0. F r (x) will have the same 
exponential tail in the left side. As expected, the Cooper pair wavefunction method leads 
to the same result for the Josephson coupling energy and the supercurrent: 

2ir 

Ej = -—< \T\ 2 > N 2 F A cos(<Pi - <t> r ) (37) 
A 

j = i-^- sin((j)i - (j) r ), (38) 



where 



* = idb (39) 



hi 



D. Shrinking of the Cooper Pair Size Due to the Insulating Barrier 

In the previous section we assumed that the tail of F t (x) in the right side, Eq. (34), is 
controlled by the same Cooper pair size, £ 0) of the left side. In other words, we disregarded 
the effect of the insulating barrier potential on the Cooper pair size. However, it is well- 
known that the Cooper pair size is reduced due to the impurity potentials. 37 ' 38 Therefore, 
it is obvious that this assumption is not appropriate especially for high T c superconductors, 
including MgB 2 , and for the insulating barriers with high tunneling resistance, R n , since 
the Cooper pair will have so much difficulty in tunneling in these cases. 

We estimate the reduction of the Cooper pair size due to the barrier potential. Equation 
(33) implies that the decay length of the Cooper pair wavefunction in the insulating barrier 
region is ~ 1/2/c. Whereas, after tunneling the Cooper pair (of size £o) will have the memory 
of the phonon-mediated attraction up to £ — d. So we suppose that the Cooper pair have 
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size 1/2k with the weight d/£ and £ with the weight (£ — d)/£o- Accordingly, the effective 
Cooper pair size £ e ff may be written as 

1 1 fr, — d d , . 

+ 2k— . (40) 



Ceff £o £o & 

Since d « (, in most cases, we find the effective Cooper pair size in the right side, 

£ ^° (41) 

Then, the effective Cooper pair wavefunction of F e (x) in the right side is 

£ _ 4 + 2/id 

F t (x) oc e ( -ff = e «o x > x a = 0. (42) 

Now, it is important to notice that this reduction of the Cooper pair size is fully meaningful 
only when the exponential factor in Eq. (42) is ~ e _1 just after tunneling, i.e., 

F,(d) oc e~«^7 = e ~ 1± ^r d ~ e -1 . (43) 



Thus, the Cooper pair size shrinks significantly when d ~ yCo/2/t, i.e., for high T c super- 
conductors, such as MgB 2 and for the insulating barriers with high R n . To put it another 
way, there is a (sample-dependent) threshold tunneling resistance, R t h, above which the 
supercurrent decreases much faster, although it is hard to determine the exact value of R t h- 
The overlap of the Cooper pair wavefunctions is given by 

F;(x)F l (x)dx = V<\T\ 2 >J2^>kJ:«^f- / e ~~ Xdx 

.V^^Eu'A^^e-^. (44) 

Consequently, for d ~ \J^ /2k the Josephson coupling energy and the supercurrent are 
written as 

27T r, r, 1 l + 2Kdj 

Ej = -—< \T\ 2 > N 2 F A——e~— d cosifr - <f> r ). (45) 
A 1 + ZKd 

and 
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1 A 1 

J = X^R n l^d e ~ 5 ° sin ^~^ (46) 

with 

1 A 1 -w d 

J1 = A^IT^ e ° • (47) 

Since the DC Josephson current is very small in this case, fabrication of MgB 2 SIS junctions 
requires very accurate nanometer scale manipulation of the insulating barriers. For instance, 
if £o ~ 100A, we need an insulating barrier of thickness d smaller than ~ 10A (i.e., R t h ~ 
0.01 — 0.1Q) to see the Josephson supercurrent, which explains why the previous experiments 
couldn't see the DC supercurrent corresponding to the big gap of MgB 2 . 14 ' 15 



E. R n Dependence of the DC Supercurrent 



We have found the maximum DC supercurrent: 

ji — T~ET for low R n (48) 



XeR 
1 A 

^ X eR, 



d for high R n (d ~ ^o/2/t) (49) 



where 



R-i = ^N 2 F < |T| 2 > . (50) 
For low R n the DC supercurrent is inversely proportional to R n , which may be called linear 



region, whereas for high R n , corresponding to d ~ y£ /2K, the supercurrent decreases more 
quickly with increasing R n , which may be called steep region, although the boundary is 
not that sharp. This behavior agrees with the experiments. 15-19 (Presumably in the steep 
region the supercurrent may drop to zero exponentially. This is the reason why we kept 
the exponential factor in the previous section.) One needs an interpolation formula that 
smoothly connects both limits, for low R n and for high R n . (Nevertheless, high R n formula 



Eq. (49) is valid in a rather narrow range where d ~ J^q/2k and the exponential factor 
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~ e x .) For that purpose, we may multiply 2nd in Eq. (49) by the approximate function 
for the Heaviside unit step function, 39 i.e., 

S n (K,d) — -[1 + tank n{2nd — 2K,d c )] (51) 

where 2nd c may be chosen to give the exponential factor equal to e -1 , while n is selected to 
have the resistance spread of factor of 10, according to the experimental data. (In this case 
2nd c gives much higher resistance than the threshold resistance Rth-) 

In the presence of ordinary impurities, the Cooper pair size also decreases from £ to £o 5 
defined by 21 ' 31 ' 38 

hhr < 52 > 

where £ is the mean free path. (It seems that Eq. (52) works better than the other expression, 
i.e., v^o- 37 ) The effective Cooper pair size is then replaced by 

ieff ~ (l + 2^(f + eo)' (53) 
It is remarkable that ordinary impurities also decrease the supercurrent by reducing the 
Cooper pair size. 

Figure 2 shows the comparison of our theoretical calculations with the experimental 
data for Pb-PbO x -Pb (PPP) at 4.2K by Schwidtal and Finnegan, 16 and Sn-SnO-Pb (SSP) 
at 1.4K by Tinkham's group (Danchi et. al.). 19 The coherence lengths £ fc> r Pb and Sn 
are ~ 820A and ~ 1800A, respectively. 40 Note that above R n ~ 40fi the supercurrent of 
Sn-SnO-Pb junction becomes larger than that of Pb-PbO x -Pb junction, due to the bigger 
Cooper pair size of Sn. For linear region ji = l.5/R n (mA) for PPP and ji = 1.2/ R n (mA) 
for SSP. McMillan and Rowell 41 observed that for a typical junction with a resistance of 30f2, 
the exponential factor is roughly exp(— 20). So the tunneling resistance is assumed to be 
R n = 30xexp(2Kd-20)(tt). The threshold resistance for PPP is R th ~ 30fi, (i.e., d th ~ 20 A), 
whereas for SSP, there are two threshold resistances for Sn and PB, leading to (average) 
value, R th ~ 5011 For PPP junction we used Sn(Kd) = 1/2[1 + tank 1.9(«d - 21.25)], 
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whereas for SSP junction we used Sn(nd) = 1/2[1 + tank 1.0(nd — 21.25)] for Pb and 
Sn(nd) = 1/2[1 + tank l.O^nd — 22.75)] for Sn, respectively. Unfortunately, the mean free 
path £ is not available for those experimental data. However, since the thickness of the film is 
about 2, OOOl, 16 we suppose £ ~ lOOOl. If we assume £ = lOOOl for Pb-PbO x -Pb junction, 
we obtain the impurity-limited coherence length £ = ££o/(£ + £o) — 450. 6 A, leading to 
d c ~ 21.23A, in agreement with the above R n (Kd). For SSP junction, the coherence length 
of Sn, ~ 1800A requires £ ~ 800A to get d c = 23. 5A, which is also consistent with the R n 
used. Considering the uncertainty in £, the perfect fit with the experimental data is rather 
fortuitous. Nevertheless, our approach explains clearly why the Cooper pair of Sn can tunnel 
through the barrier more easily than that of Pb. 

III. JOSEPHSON EFFECT AT FINITE TEMPERATURE 

Now we consider the Josephson coupling energy and the supercurrent at finite tempera- 
tures for low tunneling resistance R n . The Bogoliubov-Valatin transformations are defined 
by 42 

4t = U W] + v' k et kl 

c r H = u r k e kl - v r k e^ (54) 

: ql = U qJll ~ V qfqV 



4t = u % + yl qf- 

4 = U %1 ~ KfqV ( 55 ) 



Consequently, the coupling energy is written 

Ej = -jJ2i<v r k *(l - 2f k )u l q v l q {\ - 2f q ) + - 2/ fe )«(l - 2/,)] 

kq 

^ 2E r k x 2E\ v J A Jq> v ; 

At T = OK this expression leads to Eq. (13). After summations over k and q, we obtain 

E J = -2J ^ T ^ T \ os^ l - (t)r ). (57) 
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Notice that the pair-scattering matrix element, J, has rather weak temperature dependence 
through the effective Cooper pair sizes £ e //,i(T) and (, e ff,2(T): 

' = * < Pi' > i^e^m + * < m 2 > if^wn m 

It is well-known that the pair-correlation amplitude, or the Cooper pair wavefunction falls 
exponentially: 21 ' 26 

F(r) oc exp{ — T —) T = OK (59) 

F(r) oc exp(—£-) T -> T c (60) 

leading to the decrease of the effective Cooper pair size £ e //(T) from 7r£ /2 at T = OfT to 
3.5£o/4 near T c . (In the previous sections we used £o instead of 7r£o/2 at T = OK.) The 
limiting behavior of the effective Cooper pair size can be shown to be 43 

£e//(T)-^(l-^) near T = OK (61) 

W/(T)4^(1-^) near T c , (62) 

where £ e //(T) = 2 ^ff° A2 at any temperature. Unlike A(T), near T = OX, f e/ /(T) 
decreases with increasing temperature in a parabolic manner. 

Accordingly, for a symmetric junction the Josephson coupling energy and the supercur- 
rent are written as 

Ej = -2\T\ 2 ^-^p-Nl& 2 {T)cos<t> (63) 
TtVp A 

= l^A^ 

A ft-up eit n 



with 



A fri^ ei? n 
The reduced dc Josephson current is, then, given by 

h(T) _(u„(T)tf(T) 



A = i%S^£Q. (65) 



2 



ji(O-ft') |6 A 
A 2 (T) 



A 0V /vr 2 T 2 + A 2 (T)' 
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(66) 
(67) 



For an asymmetric junction, one finds the coupling energy and the supercurrent: 
Ej = -7T < |T| 2 > iV F , 1 iV F , 2 A 1 (T)A 2 (T)( 7r ^ l(T) + ) cos (0, - r ), (68) 

and 

J = < |T| 2 > 7V F)1 7V F!2Al (T)A 2 (T)(4f^ + ^f^) S m(0, - r ), (69) 

with 

J1 = l^A 1 (T)A 2 (T)(%iP! + %^). (70) 
2eit n Xi?ivf,i A 2 ftfF,2 

Therefore, the reduced DC Josephson current for an asymmetrical junction is 



3l (T) A X (T)A 2 (T) ^ 



Ai A2 



(71) 



3l (0K) A 1 (0)A 2 (0) f + 
Figure 3 shows the temperature dependence of the maximum DC supercurrent ji(T) 
for a Sn/Sn and a Pb/Sn junction compared with our theoretical calculation from Eqs. 
(65) and (69). Data are from Fiske. 44 As in Figure 2, we used i = 800 A and £ = 820l 
for Pb and £ = 1800A for Sn, respectively. The approximate temperature dependence of 
£ e ff(T) is given by £ e //(T) = ^/=§0+£2 > which is almost identical to the accurate numerical 
calculation of £ e //(T). 43 As can be seen, the agreement between theory and experiment is 
fairly good. 



IV. DISCUSSION 



It is clear that more study is needed to understand the intriguing properties of the 
Josephson effect. In particular, the sign of pair-quasi particle interference term, generaliza- 
tion of this approach to weak links and SNS junctions, and explanation of many unsolved 
problems in AC Josephson effect may be interesting problems. For AC Josephson effect an 
investigation based on this approach will be published separately. 45 

The Josephson effect in MgB-i requires more careful study. It is highly desirable to 
fabricate clear-cut SIS junction with very low tunneling resistance {R n ~ 0.01 — O.lfi) for 
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the detection of the Josephson current for the big gap. The insulating layer should have 
thickness not larger than 10A and a small band gap. 

Very recently, a similar problem has been found in the flux quantization of supercon- 
ducting cylinders. 46 Whereas most previous approaches focused on the phase of the effective 
Cooper pair wavefunction in the presence of the magnetic flux, flux quantization turned out 
to be due to the flux dependence of the pairing energy. 46 

V. CONCLUSION 

We have introduced Cooper pair wavefunction approach to the Josephson effect in SIS 
tunnel junctions. We have found that the Josephson tunneling depends on the size of 
the Cooper pair and its shrinking during the tunneling. Accordingly there is a material- 
dependent threshold of tunneling resistance above which the DC Josephson current decreases 
much faster with increasing the tunneling resistance. High T c superconductors, including 
MgB 2 , have tightly-bound Cooper pairs with small size, which can not tunnel through 
the insulating barrier easily, leading to extremely small critical current, in agreement with 
experimental findings. This understanding also explains why the observed DC supercurrent 
of low T c superconductors, such as Sn and Pb, decreases much faster than 1/R n above the 
tunneling resistances R n in excess of several ohms. It is of interest that ordinary impurities 
limit the supercurrent by reducing the Cooper pair size, too. We have also shown that 
the (weak) temperature dependence of the Cooper pair size contribute to the temperature 
dependence of DC supercurrent. 
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FIGURES 

FIG. 1. Cooper pair wavefunctions near the insulating barrier. The solid lines denote the 
Cooper pair wavefunctions without the barrier, while the thick lines show the change of the Cooper 
pair wavefunctions due to the barrier. 

FIG. 2. The maximum DC supercurrent vs the tunneling resistance R n for Pb-PbO x -Pb junc- 
tion (at 4.2K) and for Sn-SnO-Pb junction (at 1.4K). Data are from Schwidtal and Finnegan, Ref. 
16 and Danchi et. al., Ref. 19. Notice the crossing of the supercurrents near R n ~ 40S1. The solid 
lines are our theoretical calculation from Eqs. (48) and (49). 

FIG. 3. Temperature dependence of the DC supercurrent of a Sn/Sn and a Pb/Sn junction in 
comparison with our theoretical calculation. Data are from Fiske, Ref. 44. 
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